STRICTLY CHAINED (p,g)-ARY PARTITIONS 



LAURENT IMBERT AND FABRICE PHILIPPE 

Abstract. We consider a special type of integer partitions in which 
the parts of the form p a q b , for some relatively prime integers p and q, 
are restricted by divisibility conditions. We investigate the problems of 
generating and encoding those partitions and give some estimates for 
several partition functions. 



1. Introduction 

A partition of an integer n is a non-increasing sequence of positive integers, 
called parts, summing up to n, possibly subject to one or more constraints. 
For instance, one may want the parts to be distinct, odd, prime, powers 
of some integer, etc. The most famous reference for integer partitions is 
doubtless the textbook by Andrews pp. 

Strictly chained partitions are finite sequences of integers that decrease 
for the divisibility order. In other words, partitions of the form n = a\ + 0,2 + 
• • • + ojfc into distinct positive integers a\, . . . ,a& such that afc|afc_i | ■ ■ ■ |oi. 
Binary and more generally m-ary partitions with distinct parts are obviously 
special instances of this type of partitions. In general, an integer admits 
several such partitions; for example 

873 = 512 + 256 + 64 + 32 + 8 + 1, 
873 = 720 + 120 + 24 + 6 + 2 + 1, 
873 = 696 + 174 + 3. 

These unconventional partitions have been considered by Erdos and Loxton 
in [5]. If p(n) denotes the number of strictly chained partitions of n and 
Pi(n) the number of partitions of this type whose smallest part is 1, they 
show that p(n) > log 2 n for n > 6 and pi(n) > \ log 2 n for n > 27 except 
when n — 1 is prime, in which case p\{n) = 1. For x sufficiently large, 
they also prove that the sum function P(x) = Ei<n<iP( n ) behaves like 
cx p , where c is an unknown constant and p ~ 1.72865 is the unique root of 
Q{p) = 2, where £ is the Riemann zeta function. 
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In this paper, we investigate a special case of the type of partitions defined 
above. We consider strictly chained (p, g)-ary partitions, i.e., partitions with 
distinct parts of the form p a q b , with the further constraint that each part 
is a multiple of the following one. For sake of simplicity, we further assume 
that p and q are relatively prime integers greater than 1. 

This work arose from recent developments on so-called double-base chains, 
in particular their use in speeding-up exponentiation and elliptic curve scalar 
multiplication [3] . The first aim of this paper is to provide some theoretical 
results as a basis of further algorithmic studies. We investigate several prob- 
lems including those of generating, encoding and counting those partitions. 
We also give some results on the length (the number of parts) of the short- 
est partitions of that type. The special case min(p, q) = 2 will be given a 
special attention, in particular the case (p, q) = (2,3). In addition, many of 
our results generalize almost immediately to any number of relatively prime 
integers instead of (p, q) only. 

In the following, the set of all strictly chained (p, g)-ary partitions whose 
parts sum up to U is denoted by 0([7), and its subset of partitions of U 
with no part 1 by Q*(U). Cardinalities of these sets are denoted by W and 
W* respectively. 

2. Generating partitions in Q(U) 

2.1. Complete generation. We define three mappings from subsets of 
£l(U) to the set P of all (general) partitions; i.e., the non-increasing se- 
quences of positive integers whose sum is finite. Let ro£fl. The first map 
consists in multiplying each part of w by p, and the second one does the 
same with q instead of p. The resulting sets of partitions are denoted p f2 
and q fl respectively. 

Definition of the third mapping depends on (p, q). If min(p, q) = 2, we 
define the binary amount of a partition as the sum of all its binary parts or 
if none. The third mapping is then defined as follows: 

- if min(p, q) = 2, increase by 1 the binary amount of w, 

- if min(p, q) > 2, add a part one to w. 

In both cases, the resulting set of partitions is denoted by 1 fi. Each of these 
three mapping is clearly injective. Moreover, P 0,(U) C Vt{pU) and q Vt{U) C 
Q(qll). However, 1 Q(U) <£_ Q(U +1) in general. Indeed, if min(p, q) > 2, the 
part 1 may appear twice in 1 Q(U) (remember that our mappings take their 
values in P), and even in the binary case the divisibility condition may be 
lost. For example, the strictly chained (2, 3)-ary partition (6, 2, 1) is turned 
into (6,4) which is not in 17(10). 

If U is a positive integer the set Q(U) is possibly not empty (of course, 
if min(p, q) = 2 it contains at least the binary partition of U), whereas it 
is always empty if U £ N. We shall consider that 0(0) contains the empty 
sequence (). This convention is consistent with the simple but essential 
identities given in Lemma [T] below. The first one is readily obtained by 
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considering partitions without or with part 1, and the second one by noticing 
that, in a partition with no part 1, either p or q (or both) must divide the 
smallest part, and then all parts. Disjoint union of sets is denoted additively. 

Lemma 1. For any integer U , 

n(u) = n*(u) + x si*(jj - 1), n*(u) = P n{u/ P ) u q n(u/ q ) 

Using Lemma [TJ ground value = {(1)} and convention Q(0) = {()}, 
partitions in Q(U) can be generated recursively. 

Corollary 1. Let ko = p~ l mod q and £o = q~ l mod p. Then, for all 
U G N, 

(1) n(pqu) = p n( q u) + q (n(pU) \ p n(u)) , 

(2) n{pqu + 1) = lp n( q u) + lq (n{ P u) \ p n(u)) , 

and if 1 < r < pq 
(3) 



U(pqU + r) 



' p Q(qU + k ) + lq n( P u + p-£ ) 


ifr 


= hp 


in(pU + £ ) + l Pn{qU + q-k ) 


ifr 


= £ q 


p Q(qU + k) 


ifr 


= kp, k ^ ko 


^{qU + k) 


ifr 


= kp + 1, k / c 


q n( P u + £) 


ifr 


= £q,£^ £o 


lq n( P u + £) 


ifr 


= £q+l,£^p 


,0 


otherwise. 



Proof. Relations ^ and @ are immediate. Since p A q = 1, for 1 < k < q 
and 1 < £ < p we have W(kp) = PQ(k) and U*(£q) = q Sl{£). Moreover, 
there exist k, I such that r = kp and r — 1 = £q if and only if k = ko and 
I = p — £q- Indeed, (ko,p — £o) is the unique positive solution of the equation 
kp — £q = 1. Thus the first case in ^ is proved and the second one is 
obtained by symmetry. The other cases follow easily. □ 

These relations take the simplest form when p = 2, as the last three cases 
in ([3]) disappear. Relation ([T]) is the less efficient one in general because of 
set difference. It may be improved in this case. 

Proposition 1. If p = 2, we have for all U 6ff 

U(qU) = q fl(U) + 1 n(qU -1). 

Proof. Since g A 2 = 1, 2" ^ (mod q). Thus the binary amount of a 
partition in tt(qU — 1) is never of the form 2 n — 1, so that adding 1 to 
it does not affect the chain condition. Therefore, 1 $l(qU — 1) C £l(qU). 
Moreover, the binary amount of the result is not 0, thus l £l(qU — 1) C 
£l(qU)\ q Vt{U) . Finally, this inclusion is an equality since the binary amount 
of every partition in 0,(qU) \ q £l(U) is positive. □ 
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2.2. Encoding with words. When p = 2, relations given in Corollary [T] 
provide us with a compact representation of Q(U). An edge-labeled weakly 
binary tree with root U and labels in {1,2, q} can be used to represent the 
partitions in 0,(17) by words on the alphabet {l,2,q}. (The nodes are the 
successive arguments of f2.) Using this representation, the value of U is 
easily recovered from its encoding word: we start from the leaves (whose 
value equals 1) and go through the root U while performing operations +1, 
x 2 and x q according to the edge labels. As an example, the tree representing 
£1(19) in the case (p, q) = (2, 3) is given in Figure [T] Note that the language 
obtained this way to represent £l(U) is a hypercode (no word is a sub-word 
of another), so that other compact representations are available (the tree in 
Fig. [I] represents f2(19) as a prefix code). 




1 



6 Cj 

5 12 2 — 1 — 1 

Figure 1. 0(19) = {1112222,1112213,1332,131122} for 

(P,?) = (2,3) 

In the general case, relations Q and ^ prevent such simple, compact en- 
codings. Nevertheless, a natural graphic representation easily follows from 
the definition by noticing that strictly chained (p, g)-ary partitions are parti- 
tions with distinct parts in the set {p a q b , (a, b) G N 2 }, under the constraint 
that couples (a, b) of exponents form a chain in N 2 endowed with the usual 
product order. These chains can be encoded with words in {0, 1, 2, 3}*. Con- 
sider a length-maximal increasing path V in N 2 containing a given chain C, 
beginning at point (0,0) and ending at the maximal point in C. Such a path 
is called C -filling. Then, starting from the empty word, iteratively form a 
word by using the following rule when progressing on V: at point (a,b), the 
added letter is 

if (o, b) $ C and (a + 1, b) G V, 

2 if (a, 6) C and (a,b + 1) G V, 

1 if (o, b) eC and (a + 1, b) 6 V, 

3 if (a, b) G C and (a, b + 1) 6? or (a, b) maximal in C. 

Given a chain C in N 2 , there are generally many C-filling paths. For exam- 
ple, the chain C = (0,0), (1, 1) admits two C-filling paths, encoded by the 
words 123 and 303. In order to get an unambiguous correspondence between 
chains and words, we have to select a particular C-filling path. We choose 
the one with minimal a- values: for each i, the i th point (oj, i + 1 — dj) in this 
path is such that o« is the minimal abscissa of the i th points of all C-filling 
paths. In other words, we always go North before going East as illustrated 
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in Figure [2| We denote by wq the associated word. When C runs through 
all finite increasing chains in N 2 , it can be shown that the words wc are 
precisely those ending by 3 and with neither 02 nor 12 as factors. Note that 
the language representing (C/) is clearly an infix code (no word is a factor 
of another). 

3 t 36 36 36 



2" 



2° 



2° 



2° 



(a) 3203 



(b) 3013 



(c) 1133 



(d) 11003 



FIGURE 2. Sl(19) for (p, q) = (2,3) and encoding words 



2.3. Random generation. The tree-based representation of £l(U) given in 
the previous subsection allows for a straightforward uniform unbiased sam- 
pling method, as soon as the weight of each subtree is known. We show how 
to efficiently compute these weights in the next section. Accordingly, defin- 
ing a Markov chain is seemingly superfluous for uniform random generation. 

Nevertheless, availability of a symmetric, connected transition graph on 
£l(U) may be useful for algorithmic purposes, in particular when search- 
ing to minimize the number of parts through recoding or while performing 
addition. We construct such a graph in the special case (p, q) = (2,3). 
Transitions are based upon elementary identities which are represented in 
Figure [3] below. The second identity is a generalization of 4 = 3 + 1. 

L . L L .,, L 



2 a 2 a ~ L 2 a 2 C 

(a) 2 + 1 = 3 (b) 2(2™ - 1 + 2 n+1 ) = 3(2 n+1 - 1) + 1 

Figure 3. Valid symmetric transitions between partitions 
in Q(U) for (p, q) = (2, 3) 

More formally, choose a partition w in Q{U). For some fixed value b, con- 
sider the part 2 a 3 & in w with maximal a- value. If the part 2 a_1 3 6 also exists 



in w (Figure 3(a)), subtracting both parts and adding the part 2 a_1 3 6+1 
results in an other partition in Q(U). Indeed, the other parts 2 c 3 d in w 
satisfy either c> a, d > b or c < a — 1, d < b, so that the chain constraint 
is respected. The bounds for authorized (c, d)-values are represented with 
bold lines in the lower-left and upper-right corners in Figure [3j If the part 
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1333 



2133 



2213 



2223 



11013 13003 21003 



Figure 4. The transition graph G(27) for (p, q) = (2,3) 



2" 1 3 b does not exist in w, let C be the maximal set of contiguous parts 
of the form 2 a ~ ? 3 fe , with i running from 2 to an eventual maximal value n 



(Figure 3(b) ). Next consider the value 2 c 3 b with c = a — 2 if C is empty or 
c = a — n — 1 otherwise. This value is not a part by definition. If c > and 
if there is no part in w of the form 2 c+1 3 d with d < b, we obtain another 
partition in f2(£7) by multiplying every part in C by 3, subtracting 2 a 3 b 
from w, and adding the two parts 2 c 3 fc , 2 c 3 fc+1 . It is easy to see that the two 
families of transitions defined above are reversible, and the proof is omitted. 

Let G(U) be the symmetric graph on 0,(17) corresponding to the above 
transitions. The transition graph on 17(27) = {11013, 13003, 1333, 21003, 
2133, 2213, 2223} is given as an example in Figure |4j 

Theorem 1. For (p,q) = (2,3), the graph G(U) is connected, and its di- 
log 2 U 

ameter is at most 



log 2 log 3 



Proof. It suffices to show that there is a path between any partition in 
Q(U) and the binary partition of U. Let w be a partition in £l(U) different 
from the binary one, let b be the largest power of 3 amongst all parts in 
vj, and let a be such that 2 a 3 fe is the smallest part amongst those; i.e. 
a = min{z ; 2 l 3 6 a part in w}. Since there is no part in w of the form 2 J 3 rf 
with d > b (b is maximal), both downwards transitions are allowed. If 2 a 3 &_1 
is not in w, we apply transition 3 = 1 + 2 as illustrated below. 



2 a 2 a 



Otherwise, transition 1 + 3 = 4 (or its generalization) is allowed, as shown 
on the following figures. 
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• (•) 




(•) 



2° 



In both cases, the number of parts with power of 3 equal to b is reduced by 
at least one unit. Iterating the process a finite number of times results in 
a partition in Q(U) where the largest power of 3 is less than b, and further 
iterations eventually drop this exponent to 0; meaning that we have reached 
the binary partition. 

Regarding the diameter of the graph, it suffices to note that the exponents 
a,b > of the parts of any partition in Q(U) are located below the line 
a log 2 + 6 log 3 = log U, so that the number of transitions between w and 
the binary partition in the above algorithm cannot exceed the number of 

possible (a, 6)-values (with a,b> 0), that is, 2 io° g 2 w 3 • D 



Finding a partition with a small number of parts, in a reasonable amount 
of time, is of major importance for the applications mentioned in the in- 
troduction. For example, in the domain of elliptic curve cryptography, the 
complexity of a scalar multiplication algorithm (the elliptic curve equivalent 
of an exponentiation) based on double-base chains [3] heavily depends on 
the number of parts in a (2, 3)-ary partition of the exponent. 

If \w\ denote the number of parts of a partition w, we define cr(U) = 
mm roeO(;7) the length of the shortest partitions in £1(JJ). For example, 
from figure [2] it is easy to see that cj(19) = 2. The values of a can be easily 
computed using Corollary [T] by noticing that 

a(pqU) = mm(a(qU),a(pU)), a(pqU + 1) = 1 + a(pqU). 

The relations in ^ can be adapted as well for numbers of the form pqll + r 
for 1 < r < pq. For (p, q) = (2, 3) the following Maple code can be used to 
compute the first 500000 values of a in approximately 1 second, 
s := proc(U) 

option remember; local r; 
if U <= 2 then 1 else 
r : = irem(U, 6) ; 

if r=0 then min(s(U/3), s(U/2)) 
elif r=l then 1 + s(U-l) 
elif r=2 then s(U/2) 

elif r=3 then min(s(U/3), l+s((U-l)/2)) 
elif r=4 then min(s(U/2), l+s((U-l)/3)) 
elif r=5 then 1 + s((U-l)/2) 



3. Shortest partitions 



fi 



fi 



end: 



8 



LAURENT IMBERT AND FABRICE PHILIPPE 



As we shall see in Section 5.2 if (p,q) = (2,3) there is only one partition 
for numbers of the form 2 a 3 — 1, so that cr(2 a 3 — 1) = a + 1, whereas 
<r(2 a 3) = 1. Our numerical experiments suggest that cr(U) ~ (log 2 U)/4 on 
average, which empirically confirms an intuition that cr(U) is unfortunately 
not in o(logU), even in a loose sense. 

It is perhaps possible to significantly reduce the minimum number of parts 
by considering signed strictly chained partitions. For instance, if one allows 
the largest part in any w £ Q(U) to be less than or equal to U + 1, while 
allowing the other parts to be either added or subtracted, what is the length 
of the shortest such partition(s) of U? Studying similar relaxed definitions 
of chained partitions for the same purpose shall be the object of further 
research. 



4. Computing W(U) 

In this section we provide different formulas to compute W(U), the car- 
dinality of £l(U), with a special attention given to the case p = 2. 



4.1. Simple identities. From Lemma[T]it follows immediately that, for all 
U £ N, we have 



(4) 
(5) 



W(U) = W*(U) + W*(U - 1), 
W*(U) = W(U/p) + W{U/q) - W{U/pq). 



Corollary [T] and Proposition [T] are easily translated in the same way. Simply 
note that the sets involved in the first two cases of Q are disjoint. Let 
ko = p _1 mod q and £q = q -1 mod p. Then, for all U £ N, we have 



(6) 



W(pqU) = W(pqU + 1) = W{pU) + W{qU) - W(U), 



and, for 1 < r < pq, 
(7) 

W(qU + ko) + W(pU + p 



W(pqU + r) 



if r = kop 



W(pU + io) + W(qU + q-k Q ) 
W{qU + k) 

W{pU + i) 



if r 



if r = kp, k ^ ko 
or r = kp + 1 , k ^ 



ko 



if r = 
or r = 

otherwise. 



+ l,£^p-£o 



The case p = 2 allows for further simplifications. 
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Proposition 2. If p = 2, then for all U E N we have 
(8) W{qU) = W{U) + W(qU -I), 



(9) W{qU+l) 



inn + ir( V -1 ) /or ' '" " p - 



W(17) + W 



qU+1 



for odd U, 



(10) W(qU + r)=W 



qU + r 
2 



i/2 < r < g- 1. 



Proof. Relation Q follows from Prop. [T] at once. Now, since p = 2 we have 



Iq = \ and A;q = Then, for 1 < r < 2g, relations given in ([7]) read 



(11) 

W(2gC/ + r) 



' W{qU +^) + W{2U + 1) ifr = g+l 
W{2U + l) + W(qU +2^) ifr = g 

if r = 2k, k 4 
or r = 2/t + 1, A; / 2_i 



The last relation is summarized in (|10[ ) by considering odd and even ?7. The 
second one is a particular case of (|Sj) , and the first one is identical to Q for 
odd U. Indeed, according to we have W(qU + ^) = W(2qU + q - 1) 
and W(gC/ + 2±i) = W(2gC/ + g + 2). Finally, @ for even J7 follows from @ 
and ([§. " □ 

Note that if (p, q) = (2, 3) relations ^ may be summarized as 
(12) W(3U + 1) = W(U) + W ( 3 



C/ + 1 


-0 


2 





4.2. A general relation. Let W P (U) be the number of partitions of U 
with distinct parts taken in {p n , n £ N}. Clearly Wp(J7) G {0, 1}, depending 
whether or not U can be written in base p with digits in {0, 1} only. As an 
example, Erdos and Graham conjectured [6] that the only powers of 2 such 
that W3(2 n ) = 1 are 1, 4 and 256. (This has been verified by Vardi [7] up 
to 2 x 3 20 .) Of course, since the binary expansion of U is unique, we have 
W 2 (U) = 1 for all U > 0. It is easy to see that W p {kp + 1) = W p (kp) = 
Wp(k) and W p (kp + r) = if r mod p {0, 1}. Together with the initial 
condition W p (0) = 1, these relations define W p recursively. Those quantities 
are involved in the following theorem. 
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Theorem 2. If U N, then W(U) = 0. Otherwise, for U > 1, we have 
(13) W (U) = W p (U) + W^-j+ Yl W^)w JJ I, 



(14) 

where S Ptq (c, U) = 



i »/ LWJ 

otherwise 



1 (mod q) and Wp([/ modp c ) = 1 



Proof. Let us sort the partitions in Q(U) with respect to their p-ary amount, 
that is, the sum of all parts of the form p a . There are W(U/ q) such partitions 
with p-ary amount equal to 0. This occurs either when there is none or when 
q divides all parts. Next, note that a partition w in Q(U) with p-ary amount 
n > 1, assuming it exists, is perfectly described by its non p-ary parts, i.e., 
parts of the form p a q b with b > 0. Indeed, if c n = [log p (n)J denotes the 
largest power of p among the p-ary parts of w, each of the non p-ary parts of 
w is clearly a multiple of p Cn q, so that dividing each of them by p Cn q yields 
a characterizing partition of (U — n)/(p Cn q). Since this correspondence is 
clearly one-to-one and onto {{U — n)/{p Cn q)), the number of partitions in 
Q(U) with p-ary amount equal to n is exactly W {{U — n)/(p Cn q)). Finally, 
n is a p-ary amount if and only if W p {n) = 1. Accordingly, 



u 



W(U) = W(U/q) + ^ W p( n ) W 



n=l 



U -n 
pc„g 



Splitting the interval [1, U] into sub-intervals of the form \p c ,p c+1 ) for c > 0, 
yields 



(15) W(U) = W(U/q) W P^ W 

c>0 n=p c 



U-n 
p c q 



If p c < n < p c+1 then we have U — n> p c q only if U > (q + l)p c , so that c- 
summands vanish for c > \og p {U /(q + 1)) except for the term W P (U) W(0) = 
W P (U). 

Finally, let c < log p (U/(q + 1)). We have n £ \p c ,p c+1 ) and W p (n) = 1 if 
and only if n = p c + r with < r < p c and W p {r) = 1. Thus 



n=p L - 



C/-n 



p c -i 

£ W p (r)W 

r=0 



Moreover, p c q divides U — p c — r with < r < p c only if J7 mod p c q = p c + r. 
Therefore, 

p c -i 



^ VFp(r) VF 



r=0 



U — p c — r 



p c q 



6 p , q (c,U)W 



( 


17 




p c q 
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where 5 p , q (c,U) is equal to 1 if r = U mod p c q — p c satisfies < r < p c 
and W p (r) = 1, otherwise it vanishes. Writing U = kp c q + p c + r, it is not 
difficult to show that this definition of 5 Piq {c,U) is equivalent to the one 



given in (14). □ 



It is worthwhile to point out that if W P (U modp c ) = holds for a given 
c it also holds for all d > c. In other words, it means that the number of 



c-summands in the r.h.s. of (13) is at most equal to the the weight of the 
first digit (starting from the unit) greater than 1 in the expansion of U in 
base p. Of course, the latter remark is useless if p = 2. Nevertheless, even 
in this case many summands do vanish, as shown next. 

Proposition 3. Assume p < q and let N = \\og p {q — ■ 

For all U and c in N, if S P , q (c, U) = 1 then 8 PA {c + k,U) = for all integers 
k such that k > — c and < \k\ < N . 

Proof. Let us write 

U = (u n , . . . ,u c ,u c -i,u c - 2 , . . .,u 1 ,u ) p = [U/p c \p c + U modp c . 

If 5 Pjq (c, U) = 1, Theorem [2] ensures that there exists I such that [U/p c \ = 
tq + 1. We thus have, for each 1 < k < c, 

k 

[U/p c ~ k \ = p\tq + 1) + J> fc -V-*. 

i=l 

Moreover, each Ui with i < c is either equal to or to 1 since W P (U mod 

k P k+1 - 



p-i 



p c ) = 1. Accordingly, p + ^2 i= iP l u c -i is an integer in p 

Therefore, if we moreover assume k < N, then it is easy to see that \ U/p' 
cannot be equal to 1 modulo q, so that <5 P)( j(c — k,U) = 0. Finally, suppose 
there exists k £ [1, N] such that <5 Pi(? (c + k, U) = 1. Then, according to the 
above discussion, 5 Ptg (c, U) should then be equal to 0, which contradicts our 
hypothesis. □ 

As an example, consider (p,q) = (2,3). Here N = 1 thus at most one of 
two consecutive summands is non-zero. An extremal case is U = 4 a , because 
^2,3 (c, 4 a ) = 1 if, and only if, c is even. 

5. The sequence W 

5.1. Generalities. For any given pair p, q, the sequence W behaves rather 
irregularly, as one is easily convinced by computing its first values (see 
Fig.§. 

On the one hand, if min(p, q) > 2, then W takes infinitely often the value 
according to ([7]). In addition, it takes infinitely often the value 1 if p = 2. 
Indeed, iterating (flOj) shows that, for all integers U and a, 



(16) W (2 a q(2U + 1) - 1) = W [ qU 




Figure 5. The first 400 values of W(U) for (p,q) = (2,3) 

Therefore, for p = 2, the sequence W takes infinitely many times each value 
of the form W(qU + ^p), and in particular the value 1 = W(^-) since the 
binary partition is obviously the only (2,g)-ary partition for numbers less 
than q. 

On the other hand, we also have the following result. 

Lemma 2. The sequence W is either {0, l}-valued or unbounded. 

Proof. Assume there exists U such that W(U) > 1. Choose two chained 
partitions of U, with respective greatest parts p a q b and p a q b , and let c = 
max(o, a') and d = max(6, b'). Next construct a sequence (U n ) as follows: 
U\ = U, and U n +\ = (1 + p nc q nd )U n for positive n. Remembering the 
lattice representation of chained partitions, see Fig. [2j it is easy to show by 
induction that W(U n ) > 2 n . □ 

For the time being, we are not aware of any pair (p, q) for which W is 
{0, l}-valued. This cannot of course happen if q writes in base p with digits 
E {0, 1} only, in particular if p = 2. More generally, W is unbounded as 
soon as there exists U that writes in both bases p and q with digits £ {0, 1} 
only. But whether such an U always exists or not is a seemingly difficult 
problem, and only a positive answer would also solve ours. 

Let us next give a coarse majoration of W(U). 

Proposition 4. Let f3 6 (0,1) be the unique solution of 1/p^ + 1/q 13 = 1. 
Then W(U) < U p for U > 1. 

Proof. Let us assume p < q and show that, for all n G N*, W(U) < U 13 for 
1 < U < np. This is true for n = 1 since for 1 < U < p, W(U) £ {0, 1} 
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and U 13 > 1. Suppose this is also true for 1 < n < m and consider U € 
[mp + 1, rap + p]. Next write U = pqV + r with < r < pg. If r = 1 then 
by © W(C/) = - 1) thus W(*7) < (U - if < U 13 since U - 1 = mp. 

In the same way, if r ^ 1, first notice that the arguments of W in the r.h.s. 
of @ and are either non-integers or lie in [l,[//p], thus in [l,mp] since 
p > 2 and m > 1. It suffices then to consider the upper bounds for the 
terms in the r.h.s. of ^ and ([7]) and to check that it also holds for the l.h.s. 
term by definition of f3. For example, if r = 0, we get from ^ that 

W{U) < W(U/p) + W(U/q) < (1// + l/q P )U P = U P . 

The other cases follow easily. □ 

A direct consequence of the above is that for the particular case (p, q) = 
(2, 3) we have W(U) < U - 79 . Note that our numerical computations suggest 
an exponent close to 0.535 for the best upper bound, whereas refining the 



above method by taking into account (10) only gives approximately 0.66 



We further analyze and discuss the case p = 2 in the next section. 

5.2. The case p = 2. Although W behaves quite irregularly on the large 
scale, its local variations obey some rules, with a main pattern of length q 
when p = 2. 

Theorem 3. For all U € N, we have if p = 2 

(17) W(qU) > W{qU + l) > W(qU - 1). 
Moreover, for < r < q — 1, 

(18) W(qU + r) > W(qU + r + l). 

Proof. For 1 < i < q— 1, the binary amount of any partition in £l(qU + i) is 
clearly not 0. It may thus be decreased by 1. This transformation provides 



us with an injective mapping from Q(qU + i) to Q(qU + i— 1), so that (18) 
and the first inequality in ( |17| ) hold. 

Next consider a partition w G Q(qU— 1) with a binary amount j3 satisfying 
f3 < 2 N . As already noticed, (5 cannot be of the form 2 n — 1 with n £ N. If it 
is neither of the form 2 n — 2, it may be increased by 2 to obtain a partition 
in Q(qU + 1) because /? + 2 is still less than 2^. Note that the resulting 
partition has at least two binary parts. Finally assume that j3 = 2 n — 2 with 
n > 2. Since w 6 rt{qU — 1), we have f3 = —1 (mod q), and thus n must 
satisfy 2" = 1 (mod q). So, let k = (2 n — l)/q. Since k < 2 N ~\ w may 
be turned into a partition in Q(qU + 1) by writing k in base 2 then using 
the identity 2 n = qk + 1. Note that the resulting partition has exactly one 
binary part (the part 1). Therefore, these two operations provide us with an 
injective mapping from Q(qU — 1) to Q.{qU + 1), which completes the proof 



of the second inequality in (17). □ 



Accordingly, jumps of the function maxVF : x \— > max[/< x W(U) occur 
for certain x £ qN. We conjecture that such jumps only occur for values 
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x E g(2N + 1). For example, if (p, q) = (2, 3) the first values of max W are: 
2, 4, 5, 7, 10, 13, 17, 19, 21, 22, 25. They occur at U = 3, 9, 21, 27, 57, 81, 
165, 171, 243, 333, 345; all of the form 6k + 3. 

Nevertheless, the best we are able to prove is that, if exceptions do exist, 
they are of the form x = 2q 2 U. Indeed, first note that values x E 2gN are 
all of the form 2q(qU + r) = 2q 2 U + 2qr for some < r < q — 1. Then, 
by (§ and @ we have 



(19) W(2g 2 [/ + 2gr) = W(2qU + 2r) + PF(g 2 t7 + gr - 1) 

(20) PF(2g 2 [/ + 2gr - g) = W(2gJ7 + 2r - 1) + W(q 2 U + qr 



g + 1. 
2 



Using Theorem 3j we know that W(q 2 U + gr — 1) < W(q 2 U + gr — ^y^ 1 ) for 
all r, and also W(2qU + 2r) < W(2qU + 2r — 1) for 1 < r < g — 1 since 
g is odd. Therefore, 2q 2 U + 2gr cannot be a jump point of max W for the 
latter values of r. Hence, the only possible exception may only occur for 
r = 0, i.e., for numbers of the form 2q 2 U. We have computed the values of 
max W(U) for U < 10 6 and g = 3, 5, 7, 9, 11, 13, 15 and have not found any 
such exception. 

We have seen that W(U) takes infinitely often the value 1 when p = 2. 
In fact, we can be a little more precise, using a simple method illustrated 
below in the case q = 3. 

Proposition 5. Let (p,q) = (2,3). We have W(U) = 1 if and only if, 
either U € {0, 1} or U = 2 a 3 - 1 for some a £ N. Also, we have W(U) = 2 
if and only if, either U <G {3,4,6, 7} or U = 2 a 9 - 1 or U = 2 a 15 - 1 for 
some a £ N. 

Proof. According to @ and if U -1 (mod 3) then W(U) > 2 unless 
U is either or 1. If U = -1 (mod 3) then by ( [l6| ), we have W(U) = 
W(2 a 3(2V + 1) - 1) = W(3V + 1) = 1 iff V = 0. The same reasoning 
shows that solving W(U) = n for U only requires solving it for U = 0, 1 
(mod 3). For n = 2, then from @ and since W(0) = 1 we have W(3U) = 2 
iff W(I7) = W(3U - 1) = 1, that is, either 17 = 1 or J7 and 3U - 1 are both 
of the form 2 a 3 — 1. The latter condition implies U = 2, so that 3U E {3, 6}. 
Next, from @, we have W{3U + 1) = 2 iff W(U) = W(3[^\ - 1) = 1, 
that is, either U G {1, 2} or {/ and SL^J - 1 are both of the form 2 a 3 - 1 
with a > 0. The latter condition is impossible, thus 3U + 1 E {4, 7}. Finally, 
replacing U E {1,2} in ( |16| ), leads to the numbers of the form 2 a 9 — 1 and 
2 a 15-l. □ 



Our numerical experiments suggest that, more generally, all values in N 
are taken infinitely many times by W, but here again we did not succeed in 
either proving or disproving this. 
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6. ASYMPTOTICS 

As noticed previously, it seems rather difficult to give an estimate for the 
asymptotical behavior of maxVF. Nevertheless, it is possible to precisely 
describe the asymptotical behavior of the average of W. In this regard, our 
results are similar to those of [5|, except the fact that we obtain a more 
explicit constant in Theorem [4] below. 

Lemma 3. Let S(x) = Yli<u<[x\ W(U). Then, for all x G M+, we have 
S(x) = 2 (s ( + S (?) - S + 1 - W* ([x\). 



JPJ \QJ \PQ, 
Proof. The relation holds for x < 1 since both members are 0, with the 
usual conventions that W*(0) = 1 and that a sum vanishes if the lower 
index exceeds the upper one. Next assume x > 1. We have by Lemma [T] 
(21) 

J2 W(U) = (W*(U) + W*(U - 1)) = W*(0)-W*([x\)+2 Y W*(U) 

U=l U=l U=l 

and 

Lxj L^J 

(22) Yl W *W = E ( W ( U /P) + W(U/q) - W(U/pq)) , 
U=l U=l 

from which the relation immediately follows since LL^J/ 71 ] = L x / n J if re S 
N. □ 

Accordingly, if there exists a such that S(x) = 0{x a ) then a must satisfy 

(23) l/p a + l/q a - l/(pq) a = 1/2. 



(24) 



It is not difficult to check that (23) has a unique positive solution, which 
moreover satisfies 

a > 1 if min(j>, q) = 2, 
a = l if (p,q) = (3,4), 
a < 1 otherwise. 

As a first consequence of Theorem [4] below and ( 24 ) , the average value of 
W goes to infinity if p = 2, whereas it goes to if min(p, q) > 2 except for 
(3,4) where it goes to a constant in the order of 1. 

Besides, a is greater than j3 defined in Proposition [ZJ Indeed, (pq)^ 13 < 
1/4 since p~^ + q~^ = 1, thus 

1// + 1/q? - \/{pqf > 3/4. 



Equation ( 23 ) also reads 

(1 -p-^-^l -q- a Y l = 2, 

which is reminiscent of the equation £(p) = 2 mentioned in the introduction. 
Nevertheless, we do not try to adapt the arguments used in the quoted paper. 
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In order to cope with the recursive expression in Lemma [3j we rather exploit 
a method introduced ten years after by Erdos et al. [4] . Although the latter 
recursion does not fit the required form because of its negative coefficient, 
it may be rectified in our special case. 



Theorem 4. Let a be the positive solution of equation ( 23 ) . We have, for 
large x, 

S{x) = C m x a {l + o{\)), 
where C P) q is a computable positive constant satisfying 

-l 



In p a In q a 



P 



1 



Proof. Let f{x) = x a S(x) for positive x. By Lemma |3j we have for x > 
and i G N 



(25) 



1 



f 



p 



2 —f 



4+1 



+ 



X 



pi+lq 



+ 



1 



l-W* 



Since all terms vanish as soon as p l > x, summing both sides for i results, 
for all x > 0, in 



(26) /(.) = 1/ (5 



) + paif ( p i J + x a ^2 i 1 W 

7 i>l 1 v ^ 7 i=0 v 

Now, let fj, be the discrete measure with masses 2 q~ a at point c = log p g, 
and p~ ai at point i for each i G N*. Using (23), it is easy to check that /z is 
a probability measure, with expectation 



(27) 



E{n) 



-c + 



P 



q ^ (pa-l)2" 

Let <?(x) = f(p x )- Since / vanishes on (0,1), g(x) vanishes for x < 0. 
Therefore, and according to (26), g satisfies on M + the renewal equation 

\ - 1 

c) ■ 



(28) g(x) 
(29) 



i=Q 



p" 



i>l 

g(x — t)fj,(dt) + h(x), 

where ft(s) = £ (l - W ([^J)) 



In order to use the Key Renewal Theorem, we now check that h is 
directly Riemann-integrable (d.R.i.), that is, both h + = max(/i, 0) and 
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h~ = max(— h, 0) are d.R.i. (see, e.g., (2j p. 154]). First notice that x h-> 
Sj<x ^* ( Lp^ - *] ) i s a constant on each interval [log p n, log p (n -+- 1)), n € N*. 
Thus h , and /i + and h~ as well since p~ ax decreases with x, are continu- 
ous a.e. with respect to Lebesgue measure. Next, since W*(U) < U" with 
a > (3 > 0, 

(30) „- w < y ^ <b"J) < E _ < j^l. 

i<x P i<x y y i=0 y 

Since x \-¥ p(P~ a ) x (l — p~@)~ l is Lebesgue integrable and decreasing, it is 
d.R.i., thus h~ is d.R.i. too. Finally, h + is also d.R.i. since it is dominated 
by x i — ^ (\x\ + l)p~ ax , which is itself d.R.i. since it is Lebesgue integrable 
and decreases for x > lnp a — 1. 

According to the Key Renewal Theorem, we thus have 

lim f{x) = lim g(x) = — !— [ h(x)dx = C p>q . 

x-+oo x->oo hj{^) J R+ 

Let us next evaluate the latter Riemann-integral. 

(31) / h{x)dx< [ (N+l)^ 



rn+l 

£> + l) / 

n>0 Jn 



dx 1 

pax ~ (J _p-a)l n p"' 



Accordingly, and making use of (23), 



p a ( 2 logp <? p' 



\2 



-1 



l ' >)) C ' P ' 9 - (p a - l)lnp Q V g a + (p a -l 

I p a — 2 1 

(33) = — - r, where dip, q) = IngH In p. 

v ' ad{p,q) KF ' H ' p a H p a -l y 



Notice that, by (23) again, d(p,q) — d(q,p) = 0, so that d(p,q) may be 
replaced by the half-sum of d(p,q) and d(q,p), which yields the claimed 
bound for C Ptq . □ 
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